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Inversion  for  Mass  Spectrometer  Multicomponent  Sampling 


1.  INTRODUCTION 


This  report  addresses  the  problem  of  a  mass  spectrometer  sampling  a  neutral  gas 
consisting  of  more  than  one  component.  During  sampling,  the  component  species  are  not  onfy 
ionized,  but  also  fragmented,  thus  indroducing  mass  fragments  not  present  in  the  original  gas. 
We  wemt  to  use  the  resultant  spectral  intensity  ratios  to  find  the  density  or  partial  pressure 
ratios  of  the  original  components.  In  addition  to  the  spectral  Intensity  ratios,  the  following  are 
assumed  to  be  known:  the  fragments  produced  from  each  original  component:  the  branching 
ratios  of  each  of  these  fragments:  and  the  ratios  of  ionization  efficiencies  for  each  component. 

In  principle,  all  the  known  quantities  are  independent  of  particular  mass  spectrometer 
characteristics,  such  as  multiplier  gain  for  example,  so  that  the  procedure  need  not  take  into 
account  the  peculiarities  of  a  specific  instrument.  It  is  required,  however  that  the  branching 
ratios  and  hence  the  electron  beam  energies  be  the  same  for  both  branching  ratio  and  sample 
measurements. 

Note  the  frequent  reference  to  ratios  of  values  rather  than  values  per  se.  This  is  because 
m£my  of  these  values  vary  from  instrument  to  Instrument,  and  further,  for  any  given 
Instrument,  with  time.  As  noted,  we  are  seeking  an  inversion  that  is  Independent  of 
Instrument  properties,  such  as  ion  source  characteristics  or  multiplier  gain.  If  these  properties 
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are  assumed  to  remain  fixed,  then  Inversion  Is  quite  simple.  If  they  are  not  fixed,  then  the 
method  described  here,  a  more  complicated  procedure,  or  possibly  some  other  procedure  must 
be  resorted  to. 


2.  GENERAL  SOLUTION  OF  THE  PROBLEM 

Let  Cj  be  the  density  or  partial  pressure  of  the  neutral  compounds,  m  In  number,  j=l  to  m. 

Let  ttj  be  the  Ionization  efficiency  of  the  jth  component. 

Let  the  By  be  the  spectral  intensities  of  the  1th  fragment  of  the  jth  compound.  This  will  be 
an  array  of  n  rows  (I  runs  from  1  to  n)  and  m  columns  (j  runs  from  1  to  m).  The  rows  identify 
fragments  of  a  given  mass  (Including  In  particular,  the  mass  of  the  unfragmented  Ion):  the 
columns  Identify  a  given  Input  component. 

Let  I=I  0)  Index  the  mass  of  the  jth  unfragmented  component,  so  that  By  when  I  Is  written 
without  an  argument  Is  a  one  dimensional  array,  dependent  on  j  only.  Bi(/)j,  however  (l(/) 
written  with  an  argument)  Is  once  again  a  member  of  a  two  dimensional  array,  with  subscripts 
l(^)  and  j. 

At  this  point,  to  clarify  the  notation,  we  illustrate  by  displaying  a  specific  example.  For  this 
purpose  we  choose  a  gas  having  four  components:  CO.  CO2.  C2H4O2.  and  O2.  and  assume 
that  the  first  two  and  last  of  these,  but  not  the  third,  Is  decomposed  into  every  possible 
fragment.  Let  the  compounds  be  sequenced  in  the  order  given  above.  Then,  the  array  for  the 
By  looks  as  follows.  (M  is  the  fragment  mass  number). 


Table  1.  Illustrative  Table  of  Symbolic  Values  for  By. 


M 

Cl 

(CO) 

C2 

(CO2) 

C3 

(C2H4O2) 

C4 

(O2) 

12 

Bii 

Bi2 

Bi3 

Bi4  =  0 

16 

B21 

B22 

B23 

B24 

28 

B31 

B32 

B33 

B34  =  0 

32 

0 

II 

PCI 

B42 

B43 

B44 

44 

0 

II 

in 

B52 

B53 

B54  =  0 

60 

Bei  =  0 

B62  =  0 

B63 

B64  =  0 

In  general,  most  of  the  By  will  be  zero,  but  for  this  relatively  small  system  they  are  not.  We 
have,  m  =  4.  n  =  6.  and  1(1)  =  3:  1(2)  =  5:  1(3)  =  6:  1(4)  =  4. 

This  example  will  be  used  again  in  Section  3  to  demonstrate  a  sample  solution. 
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Returning  now  to  the  general  development,  we  have  the  component  conservation  relations 


(la) 


“2^2"  ®12 


(lb) 


(10 


(Id) 


a. 


Let  0  =  ^  and  B,  =  .  The  p.,  are  relative  branching  ratios,  and  the  p 

By  J  ttj  «  J 

efTiclency  ratios  referred  to  in  the  introduction. 

C, 

Lety,  =  — i- 


the  ionization 


i  C, 


J+i 
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Then 


^2  P2  ®12^  ®i3  P2  B^  ^  ^12'^  ^  P 

1=1  1=1  ^j3  i»l  i»l 


n  n 


B..  n 


^m-I  *  Pm-1  2  ®l.n.-l/ S  =  ^m-l  S  '’>■” 

1=1  1*1  1*1  1*1 


Let  Bij/Bij+i  =Xj  :  j=ltom-l  :  Xn,=  l 


i  pij/  i  p,.j+i = Sj- 

1-1  j  1=1  •'  ^ 


Then  yi  =  PiXjSi 


y2  =  P2X2S2 


yk  ~  Pk^kSk 


ym-l  ~  Pm-l^m-lSm-l 
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The  reason  for  the  condensation  of  notation  Is  that  all  pk  aJid  Sk  are  known.  We  now  need  to 


obtain  the  Xk  to  complete  the  solution. 


Let  Qj 


m  m 

®I(2)J 


m  m 

m  m 

“  jlj  ®I(k+l)J 


n 


in-1 


m  m 


The  are  the  intensity  ratios  of  the  MEASURED  SPECTRAL  LINES,  having  the  masses  of 

the  unfragmented  input  compounds.  We  now  show  that  the  right  hand  sides  of  the 
expressions  forOj^  cem  be  expressed  in  terms  of  the  and  Xk-  Since  theQj^  are  m-1  in 

number  these  then  constitute  a  closed  set  of  simultaneous  Implicit  equations  for  Xk.  When 
solved  for  the  Xk  they  give  us  the  solutions  for  the  required  values  of  yk. 

We  could  just  as  easily  have  selected  any  other  of  the  (m-1)!  sets  of  spectral  ratios  formable 
from  the  Intensities  of  the  m  unfragmented  species.  However,  the  choice  made  is  readily 
identifiable  and  seems  a  natural  choice.  The  manipulation  of  the  expressions  for  to  put 

them  in  terms  of  and  Xk  will  be  done  for  k  =  1 .  The  other  relations  then  follow  by 

augmentation  of  the  k  index. 


_  ^1(1)1  ^1(1)2  •  ^KDJ  •  ^i(l)m  (4a) 

^  ®I(2)1  ®1(2)2  ■*■•••  ®l(2)j  •  ®I(2)m 


^11  ^12  Pl(l)2  '*'•  ••  ^11  Pl(l)J 

Pl(2)l  ®I2  Pl(2)j 


'  ^Im  Pl(l)in 


(4b) 
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Note  In  the  subscripts  of  B,  the  switch  from  argumented  I  (1(1),  1(2))  in  Eq.  (4a)  to 
UNARGUMENTED  I  in  Eq.  (4b).  This  is  preparatory  to  expressing  ratios  of  By  In  terms  of  xj. 
Next,  divide  numerator  md  denominator  by  Bim  and  observe  that 


B,,/B,  =  x,x.  ,  .  .  •  x^ 

ly  Im  j  j+1  m- 


m-l 

1= n 

H 


(5) 


so  that 


m  m  m  m 

=  z  pi(i)j  n  £  pi(2)j  n 


(6a) 


(In  these  as  well  as  the  remaining  sums,  we  have 


Pl(l)l  “  Pl(2)2  -  •  •  •  Pl(j)j  -  •  •  •  Pl(m-l),m-l 


In  like  mrmner. 


m  mm  m 

“2=spi(2)jn  viPiouH^i 


m  m  m  m 


m  m  m  m 

£  Pl(ni-l)j  n  *i/ £  Pi(m)j  0 


(6b) 


(6c) 


(6d) 


III 

Letz  =  n  *!•  1  ^  ai. 

^  i=J 
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Then,  clearing  the  kth  equation  of  fractions 


Pl(k+l)J 


Ill 

Zj  =  I  P 

J  j-l 


l(k)J 


(7) 


Let  =  12^  Pl(k+l)j  “  Pl(k)J’ 
Then 


y  A,,  z,  =0  1 S  k  ^  m-l  (8) 

jtl 


form  a  set  of  m-l  linear  equations  for  the  m-l  unknowns  zj  which  can  be  solved  by  one  of  the 
usual  methods.  Following  this,  the  xj  are  recovered  from 


1  ^  j  ^  m- 1 . 


This  completes  the  solution. 


3.  EXAMPLE 


To  illustrate  these  Ideas,  we  first  assign  values  to  the  array  for  the  By  shown  In  Table  1  and 
also  to  the  required  values  of  aj  the  Ionization  efficiencies.  Here,  the  units  for  the  values  of  By 
(the  spectral  Intensities  of  the  1th  mass  fragment  from  the  jth  compound)  may  be  taken  as 
microamperes:  euid  those  for  the  Ionization  efficiencies,  when  normalized  to  a  pressure  of  one 
Torr,  as  microamperes  per  Torr.  Note  that  the  column  sum  of  spectral  intensities  is  set  equal 
to  the  corresponding  value  of  aj.  This  assignment  is  somewhat  arbitrary.  Strictly,  the  sum  of 
the  branching  ratios  for  each  component  is  unity,  and  dividing  each  value  of  By  by  oj  will  yield 
this  result.  If  the  By  are  the  actual  fragment  spectral  intensities  due  to  an  arbitrary  mixture  of 
gases,  neither  of  these  is  true:  each  column  will  be  multiplied  by  an  (unknown)  constant, 
dependent  on  the  partial  pressure  of  that  particular  component.  Only  with  this  assignment 
will  Eq.  (1)  be  true.  These  constant  values  however,  are  precisely  those  we  are  seeking  with 
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this  treatment,  and  so  are  not  known.  However,  since  all  equations  except  Eq.  (1)  deal  only 

with  RELATIVE  values  for  all  variables,  the  actual  values  as  well  as  the  units  turn  out  to  be  of 

no  Import.  Only  the  RELATIVE  VALUES  In  each  column  are  required.  To  the  right  of  the 

array,  are  listed  the  intensities  of  the  MEASURED  mass  spectrum.  From  these,  will  be  formed  , 

the  three  (=  m-1)  values  of  required  for  solution.  We  note  that  none  of  the  values  given  here 

are  actual  measured  values:  those  for  CO,  CO2  and  O2  are  reasonable,  but  all  have  been  ^ 

fabricated  for  purposes  of  illustrating  the  method  of  solution. 


Table  2.  Values  of  By. 


Cl 

C2 

C3 

C4 

BHB9 

(CO) 

(CO2) 

(C2H4O2) 

(O2) 

12 

0.0025 

0.0009 

0.005 

0 

0.002 

16 

0.05 

0.03 

0.01 

0.15 

0.2 

28 

0.5 

0.06 

0.02 

0 

0.4 

32 

0 

0.0009 

0.005 

1.5 

1 

44 

0 

0.3 

0.02 

0 

1 

60 

0 

0 

0.2 

0 

2 

ai 

0.5525 

0.3919 

0.260 

1.65 

Next  we  form  the  table  of  the  py  by  dividing  By/By.  We  also  form  the  .  and  the  Pj. 


Table  3.  Values  of  py,  Oj,  and  Pj. 
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From  these. 


4 


S  =  UQS.  =  0.8271; 
1  1.360 


j  =  h.360  =  1 046; 
2  1.300 


S  =  1,3QQ.  =  1.182 
3  1.100 


emd 


yi  =  0.5864  Xj; 


y2  =  0.6935  X2: 


ya  =  7.506  X3. 


To  complete  the  solution,  we  next  form  the  equations  for  the  x^. 


a  =  x,x^  +  0.2x^3  +  0.  1X3  +  0  ^  x,x^-t-0.2x^  +  0.1 

1  0  +  X2X3  +  0.1Xg  +  0 


X2+0.1 


= 


.,0.1 


0  +  0  +  Xg  +  0 


= 


O  +  O  +  J^  +  0 


0  +  O.OOOx^Xg  +  0.025Xg  +  1  O.OOOx^Xg  +  0.025Xg  +  1 


Here,  the  system  is  so  simple  we  need  not  introduce  the  zj;  instead  we  find  directly 


xi  =  0.05;  X2  =  0.40;  X3  =  2.1 1 1 


and 
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yi  =  ICOI/ICO2I  =  0.02932 


72  =  ICO2I/IC2H4O2I  =  0.2774 


♦ 


73  =  IC2H4O2I/IO2I  =  15.84 


4.  THE  DEGENERATE  CASE 


It  may  turn  out  that  two  (or  possibly  more)  of  the  unfragmented  input  compounds  have 
identical  masses.  A  good  example  of  this  would  be  the  presence  in  the  measured  sample  of 
CO2  and  N2O.  Since  this  only  appears  as  a  single  peak  in  the  measured  spectrum,  we  would 
then  seem  not  to  have  the  required  m-1  values  ofOj^  needed  to  complete  the  solution.  If 

however,  the  sum  of  the  components  is  treated  as  a  single  component,  we  can  still  obtain, 
instead  of  the  m-1,  m-2  ratios  of  the  input  components,  with  the  sum  of  the  degenerate 
components  appearing  cis  a  single  component. 


5.  UNRESOLVED  ISSUES 


We  have  tacitly  assumed  that  real,  positive  solutions  for  the  can  always  be  found.  This 
seems  to  be  intuitively  clear  on  the  following  grounds.  If  instead  of  inversion,  we  start  from  any 
known  mixture,  we  will  always  produce  a  specific  spectrum  containing  both  unfragmented  and 
fragmented  spectral  lines.  Thus,  when  the  process  IS  inverted,  we  know  that  a  solution  exists; 
it  is  reasonable  to  assume  then,  that  this  solution  will  be  obtained.  Still,  such  reasoning  does 
not  constitute  a  proof,  and  no  proof  has  yet  been  found. 

If  this  reasoning  is  correct  however,  we  have  here  a  meeins  of  checking  on  the  quality  of 
experimental  data.  If  a  solution  cannot  be  found,  then  this  indicates  that  there  is  an 
Inconsistency  between  the  spectra  measured,  the  measured  branch  ratios,  and  the  measured 
ionization  efficiencies.  One  aspect  of  this  question  might  be  pursued  further  by  Investigating 
the  way  assumed  errors  in  the  measured  spectra  are  propagated  backward  into  the  calculated 
original  components. 

Another  question  is  whether  it  is  possible  by  further  work  to  find  the  ratio  of  input 
compounds  even  when  they  do  have  the  same  mass.  This  issue  is  not  pursued  further  here. 
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